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A. Problems 


None 

B. Accomplishments 

Collection of ground data in support of HCMM overpasses and for 
calibration of the finite-difference heat flow model has begun. 

The finite-difference model was adapted for use on an Hewlett- 
Packard 9825 minicomputer (See Appendix A). 

A series of model calculations were made to simulate effects of 
near surface soil moisture (within 50 cm of soil surface) and water 
tables at deeper depths on the soil depths in soil surface temperature. 

During the next reporting period, collection of ground data will 
continue. Analyses of HCMM data will begin upon receipt of the data. 

C. Significant Results 

In order to investigate the general relationship between surface 
temperature and soil moisture profiles a series of model calculations 
are being carried out. Soil temperature profiles are being calculated 
during a complete diurnal cycle for a variety of moisture profiles. 
Preliminary results indicate the surface temperature difference between 
two sites measured at about 1400 hours is related to the difference in 
soil moisture within the diurnal damping depth (about 50 cm). The 
model shows this temperature difference to vary considerably throughout 
the diurnal cycle. 

The presence of a water table below the diurnal damping depth was 
simulated by fixing the soil temperature at 50 cm depth to be 2° C 



cooler than a similar site without water. Model calculations comparing 
two sites with and without this water were made. Preliminary results 
show the presence of such a water table results in a surface 
temperature difference which is constant throughout the diurnal cycle. 

If the two sites also differ in soil moisture in the top 50 cm layer 
of soil, the presence of a water table at one site simply alters the 
surface temperature difference by a constant amount. Thus the model 
may allow one to distinguish the effects of soil moisture near the surface 
from those due to a water table below the diurnal damping depth by 
consideration of surface temperature differences noted by HCMM at its 
day and night overpasses. 

It should be noted that these series of calculations assume the 
same heat flux into both sites and do not account for crop canopies. 

These factors may complicate use of the model in actual field conditions. 

D. Public at ions 
None at this time 

E . Recom men dati on s 
None at this time 

F. Fu nds Expend ed 
$29,007.94 

G . D ata Util i ty 

HCMM data has not been received 


APPENDIX A 


Soil -Temperature Model Program: 

This program is set up to provide a model of soil temperature 
for 50 cm of soil. This is broken into fifty (50) equally spaced 
layers. 

PARAMtTERS: 

X(l) is the distance between nodal points in centimeters. 

T is the time interval between calculations in seconds. 

N is the number of equally spaced notlal points. 

P is the time between print-outs in seconds. 

M(l) is the initial starting minute. 

H(l) is the initial starting hour. 

M(2) is the minute when the calculation is to end. 

H(2) is the hour when the calculation is to end. 

Y(1 ) is the ending day. 

W(l) is the adjusted air temperature. 

A(I) is the temperature of Site A at depth I. 

B(I) is the temperature of site R at depth I. 

0(1) is the moisture of site A at depth I. 

P(I) is the moisture of site B at depth I. 

Q is the amount of soil. 

K(0,l,2) is the conductivity of air, water, and soil respectively. 
G(l,2) is the heat capacity of water and soil respectively. 

G(0) is a dipolarization factor as defined by deVries 
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Program Listing 


0; 70r>-»r0;TtO 7<^ 
*1101''. 


* Subroutine "QHEAT" 

This subroutine is used to calculate the surface heat flux throughout the day. 


1: " n 1 r7 1 

2 : >Ul|+r2Wr?2 
1: r2 2-('.. 7S-r’n 


Conversion of Temperatures in Degrees Centigrade to Degrees Kelvin. 
The adjustment factor ( /lOO) allows for larger precision in the 
calculation. 


' : (’•'( n+''7 2. K, ) /l0':-»r?4 

: ( \f01 +27 1.1 n /10G>r2'i 

; { 'M 0 1 + ? 7 1 . 1 G ) / 1 ''' 0 -► r 2 

* 1 1 1 2 


Calculation of Surface Heat Loss Term 

(.000136 corresponds to the Stefan-Boltzman constant) 


7: . nnoi K.*r2'' “'1 ■*r'57 

.1 : - ( . 0 1 G * r " 4 - r 7 ) - r 1 1 

G: - ( .nO'M lG*t?r ‘4-r'* 7) »r I 2 
* 1 '>07 
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Calculation of Surface Meat Flux 

This test to determine whether the time is at night or during the day 
and adds a solar heat flux term for the day. 
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er detail see plottin^T subroutine packaj^e). 

nvo" ; 

Wtb rn, 27, r, int((n1-y)(,/f..n, ip.t ((^ 1 _v);,,)^i„, . 

rot I I , I ^ [ [ ' ) /' ■ • ) , J t - t ( ( : , V - V ) - J 

" i T'OUrJ '* • 


1 r> ; •' rrovo" ; 

T 7: 

1 : rot 

1 ^ r " i n'ovr " ; 
7<7 : wtb t b,? 
21 


rot 


?2 : 

" r> 1 * •• ; 


23: 

' t b t , 27 , <7S, i pt ( ( •''1 

-X) L 

2 1 : 

i' *>7 = 0 ; 4 *p3 


2*^ : 

i 'O’-IT) ; 1 (1, 27 , ' 

’ , 1 . 

2r, : 

wt‘ r(',T>7;Ktr- rn,v; 

f 9 

27 : 

if ''3-4 0;wtb ri',27,w 

‘ ! 1 

?'l : 

rot 

’ 9 ' t 

2‘» : 

"chor" : 


^ 2 : 

if n? = 0;S»n?; 0*0.7 


n : 

'. tb n', 27 , 4''.,r'l , iot {' 

/o 

12 : 

rot 


3 3 : 

" pv- }7." : 


74 ; 

o1 *1) • 


I*". : 

vi-b r0,-7,?q, int (^40 1 

1 9 It/.- 

1 '7 : 

rot 


17 : 

f^c.l " : 


T 4 • 

1 7,)i.v (nO_^^ ^ 


30 : 

OfiM / ( T^.J - f. 7 ) 


•;n ; 

o 1 * Y ; r: 7 * V 



■ 11 : r r> t 

’2: "x-xi:-'*: 

’7: vt-u rP,27,4''.,qt^,n,';,n 

4J : ir o4 -.4 . .’o ? 


ORICrlNAi: PAGE IS 
OF POOK QUAUTYi 


1 ■? 


1 o . r , 

\ t * ' r 0 , 2 7 , . i p t' f ^ — 'M 1 / • \ 

t 0, Z’ -i ; v.’t '• t 1. ^ ’')/' M . 

/. 


*T ^ t 0, /> ■> ; v.’t , r ^ 

r., -.7, n.,,i:„,., _ 


r r. t 


< 


0 0 


! . 

•■ V.' 

’ i r “ • 

1 . 

1 • 

V-; f ] 

' f0,27, ir,, 17 t ^ i 

•> . 

i ' 

” 5-- 1) r. 'vi t;*i - 't . >• 

' ^ : 

i f 

o-’r.O . ,>.l 1 ^n;i ' 

.4 ; 

\ , t '■ 

'■ <1 / z' , 0 s , >■ o t ( ( 


o7 ♦ 

’'0;wth rP,.*^;vt 

1 : 

wt !■ 

» (! *' 7 , 1 1 

' 7 : 

r o t 


' : 

" r- 

,-('P •' . 

1 : 

i f 

ol < n ;,1 t- o +.7 

0 ; 

r t >- 

ri\ 37 ; ir n 2 ( ( '• 1 

J 

V t! 

^ 0 . -j...,. n ^ 

n . 

r o ( 


■' ; 

■■ oj 

i - " . 

1 . 

i '■ • 

" 1 < 0 ; ■! r o +■? 

: 

\ *t b 

rP , ) ; ir r 7 { ( O ) 

• 

, ♦ u 

1 :; j f,- 


rot 


; 

ct 

• ' ** • 

\ . 

1 ' 

7 7 

j 

• r ' 

’ '\27, ■ t 


’/’''M( l-viw) 
* ► .) > - ■ 


- I ; f. ./ 




■ ) / ’ ■ • ( ( ■ i- . iy) 
( > t * ' ) •' 


1 ) - ■ ■ 


' ' ' ' 2 , 7 ^ t { 1 ♦ 1 / , j ^ , 

■ ' ' ' ' / 27 , 7 i-> , i n I ( J - , > 

^ . \. t ! 70 ^ in I ('V. ‘ 5 ^ 

/ T » * M ^ ’ 


} 

"1 * 


'^2 


The Main Program Begins Here. 

Parameters used are as previously defined and are called for bv the computer. 


: r’ir a I 0 ; '•'(I I , MO : '■>n 1 , ' I 0 ; «\1) I , 1. 1 1 , M 1 I ' : ' 1 

7 ; 'Mtr { 1 i , V ( 0 ; 1 0*^. ) , W M 1 

1 : f’iir ‘MO :4 ) ,M (0 :?] , K jo :2 I ,'MO : .M , 1 |0 : SO) ,'M'- : 1 1 , ' 10 : 

1: ont v|l ) ,T,N,P,Min ,'Mll ,Y(1 I ,IM?) ,'M?1 , Ml 

1; -1*I,(l) 

!: for 1*'2 to sn ; [. | I- 1 ] - 1 * i. j i ) ; 0 *^x 1 I 

n-0*‘M31*^'(41-*rlO 

t ; M n 1 *S0f'’M ' - M 

1 : 1 ♦ r 1 S 

>1^07 
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* Specification of Initial TeMperature Profile 

A(0) corresponds to the surface, H(I) 1 cm etc. 

>: out A(0) ,A(11 ,A(S) , MIO] ,‘>|:S] ,A(S71 ,^M^M Ml M M '' 1 , ' I 1 M , .[ 7 ; \ ,'M‘»l ! 
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: A r n - « n 1 * M M 
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Testing for Calculat ion of Conductivity and Specific Meats. 
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Specification of Conductivity and Heat Capacity Profiles. 
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* Call Subroutine to Calculate Surface Heat Flux. OF FOOR QlJ*VU.r*i 
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Calculation of Internal Nodal Temperatures. 


^!3: '■p;. 1 = 1 t o J- 1 ; ( [ I j / r 1 -C f 1 - 1 1 -2 * •.' ( I ’ -f' [ 1 + ^ M ‘ f i ! * r . 

: ( ( ■• ’ f I -1 H ■■' ( 1 1 ) * M I- 1 1 + f '' ( 1 1 +C [ M 1 1 1 * •' I H M ♦ i"' ) M 1 / ' , i I 1 * I I ] 

f P ! i ) /r 1 -I) I I-l 1 -?*D [ I ] -P ! I + i 1 ) * '• I 1 1 *r 7 

‘ : ( ( r. [ I -1 1 + 0 I T 1 ) *'M I-n + (HI T I +D( I + n ) * I I + ’ I + r7 ) * ( 1 / !• I 1 I * q I I ; nt.fX t 1 

’ 7 1 ] n 


Calculation of Surface Nodal Temperatures. 
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■^9 : 1 1 1 * r n 4 (C ( 1 1 /r 1-?*^ I n-?*C! ? 1) • !0 1 + 1 p ’^ : (;’| ) * \ I ! I > 

40: rl/n(l) *r(01-P{0) 

4 1: 4*X I n *r 1?4(F n 1 'rl-.?*nl 1 1-2MM. ?1 ) ‘P (ni + (T. n( 1 1 42*0(2) ) *» (1 P 
42: rl P’ 
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Calculation of Lower Boundary Temperatures. 


4 3: A I N P r ( N ) 
4 4 : 0( N 1 ♦ F ( N 1 
1 0319 


Reassignment of Nodal Temperatures for Next Iteration. 
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Testing for Print-out Time 
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Printing of Pertinent Data 
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Allocation of Three (3) Equal Columns for the Data 
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int(N/3) 
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End of Print -Out Routine 
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This Portion Plots the Data, 
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Testing if Calculations have run for desired time 
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Calculation of Difference every 20 minutes 
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This portion \k>ts the "5 cm Temperature Difference”. 
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